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O Abstract — Two general empirical formalisms are presented for the spectral density which take 

^into' account the deviations from the Lorentz line shape in the wing regions of resonance lines. 

These formalisms satisfy the detailed balance requirement. Empirical line shape functions, 
j which are essential to provide the continuum absorption at different temperatures in various 
* frequency regions for atmospheric transmission codes, can be obtained by fitting to experimen- 
; tal data. 

i 



I. INTRODUCTION 

Transmission codes 1 "* which calculate atmospheric transmittance and radiance for various ge- 
ometries and atmospheric conditions are widely used in remote sensing, communication, etc. It is 
essential for the further development and refinement of these transmission codes to take into 
account the continuum absorption by H 2 0-H 2 0 and H 2 0-N 2 pairs in the window regions of the 
Earth’s atmosphere. Although a lot of effort has been made for decades, both experimentally and 
theoretically, 7 ' 12 we have not yet achieved this goal. However, it has been shown that the results 
of line-by-line calculations using the Lorentz line shape, which is valid near line centers, cannot 
be used to explain the continuum absorption; this implies that the Lorentz line shape must be 
modified in the wings of the lines. 

Recent theoretical studies 13 " 21 on continuum absorption have led to significant progress in our 
understanding of the physical mechanisms and enable us to predict its magnitude and its 
temperature dependence. These studies complement existing empirical models especially at low 
temperatures at which adequate experimental measurements are lacking, however insufficient the 
present accuracy may be for atmospheric applications. On the other hand, empirically modeling 
the continuum absorption by introducing an adjustable parameter, a x factor, 22 to modify the 
commonly used Lorentz line shape provides a simple methodology for practical applications and 
it is widely used in transmission codes. 1 " 6 Of course, the reliability of the empirical models is limited 
by various other factors. One of them stems from uncertainties in the available experimental data 
on which the empirical models are based. Another is the limited temperature range and frequency 
regions covered by the reported measurements. In addition, the form of the empirical line shape 
correction affects its accuracy in representing the continuum absorption. Unfortunately, the 
empirical model 23 that is widely used in trans’fnission codes fails to satisfy the principle of detailed 
balance which governs the emission and absorption processes. We focus our attention on the 
detailed balance requirement in the present paper and show how the empirical model can be 
improved. . 

In Secs. 2.1 and 2.2 we briefly outline the basic formulae for the correlation functioif and the 
spectral density and give explicit expressions involving the Lorentz line shape. We analyse the 
spectral density obtained with a x factor modifying the Lorentz line shape and identify in Sec. 2.3 
the source of its failure to satisfy the detailed balance requirement. Then, in Sec. 2.4, a simple 
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improved empirical model which satisfies detailed balance is presented. We discuss in Sec. 2.5 a 
more general empirical model in which not only the detailed balance principle, but also the 
asymmetric character of line shape has been taken into account. Finally, some conclusions are 
drawn in Sec. 3. 


2. THE GENERAL FORMALISM 
2. 1. The correlation function and spectral density 

As is well known, the absorption of radiation at frequency ct>(cm~') by a gaseous sample in 
thermal equilibrium at temperature T is characterized by the absorption coefficient ct(a>) 

«(«) = Wc a ' thiha}l2kT) h {* e-n</* (H, (0) • fi (H K »> + <rt) • /i (H, (0)>] d t. (I) 

In this expression, fi (H \t) is the total dipole operator of the gas in the Heisenberg representation 

( 2 ) 

and the angular brackets denote the ensemble average (i.e., the correlation function) 

< /I (H)(0) . ,,<«>(,)> . £ p,< i | p m \0) • fi m, (t) | «>, (3) 

l 

where p, is the normalized Boltzmann factor of the initial state i. The spectral density F(w) is a 
Fourier transform of the correlation function defined by 

F(oi) = -!* e- to 'O«(0) • ft(t)> dr, (4) 

and F(—a>) is defined by 

F( ~ <0)s 'k} *~'"<i‘(*)-m>dt. (5) 

For simplicity in the above expressions we have dropped the Heisenberg superscript. The spectral 
density must satisfy the detailed balance requirement 

F( — cu) = e~ Mtr F(u}). (6) 

We note that the definition of the Fourier transform in Eqs. (4) and (5) is the same as that used 
to derive the empirical model 22 - 23 but differs from that used in our previous papers 17 * 21 by a minus 
sign in the exponential. 


2.2. The correlation function and the spectral density corresponding to the Lorentz line shape 
Within the impact approximation the correlation functions are given by 

<R(0) • M/)> = I Pi\d fi \ 2 e-* w+ “>\ (7) 

if 

<M0 ■ A*(0)> = I Pi - O«(0) • /.(/)>*, (8) 


where a> fi = -cu (/ , y fi = y i} and \p n \ 2 = |/i, r | 2 , and the spectral density is given by 




F(q)) = — £p*Im/I 2 

27C 


=-l p,m 2 




{u -%) 2 + y)’ 


(9) 
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and 




1 


Vfi 


Jt + cu*) 2 + y} 


1 


Vfi 


( 10 ) 


= e -fc u /trly e *« u - lu> )/*r , ,2 

*7 ^ (cu-cu^ + yjT 

The fact that the line shapes obtained based on the impact approximation are Lorentzian is clearly 
shown in Eqs. (9) and (10). The spectral density given above does not satisfy the detailed balance 
requirement. However, in the band center region where most of the contributions come from 
nearby strong lines and where the Lorentz line shape is valid, this requirement is approximately 
satisfied. The symmetric spectral density can be expressed as 


F(cu) + F(-cu) = ^i X 


r_ 


Vfi 


Vfi 


4?r ^ 0 /2£ T) (_(oj -cu^ + yjs (cu + a)*) : 

-Hal y § F y * i y * 1 

4n 2 n a ~ 0 l(cu-cu*) 2 + yj (cu +£D / s) 2 + yjJ' 


+ vi 


where the usual line strength S fi is defined by 

47T^ 

and a reduced line strength is defined by 


( 11 ) 


( 12 ) 


&■ 


0) fi th(haj fi l2kT) 


(13) 


2.3. A widely used empirical model of the spectral density for the continuum absorption 

It is well known that the Lorentz line shape is not valid in the line-wing region. In order to take 
into account the derivations from the Lorentz line shape, an empirical x(cu) factor has been 
introduced 23 in the expression of the spectral density, Eqs. (9) and (10), and it has been arbitrarily 
assumed to be an even function, 


^(“) = -2>,m 2 p V 2 x(o>-oi fi ), (14) 

* t (co-c o fi r + y} 

F(-cu) = - X PiM 2 , , fy; — x(a> + <%)• (15) 

nf (co + c o fi y + y} 

Then, the symmetric spectral density is given by 


F(cu) + F(— cu) = I — 7l 7J~ nbf \ 
4n « J, > o ‘ °fi th (ho)fi /2kT) 


Vfi 


(cu + c of) 2 + yj 


X(w + c Of) 


r — 

(_(CU -CU*) 
] 




3 he 1 _ . f ■■ y# 

4n 2 n Jo * .(w - cu*) 2 + yj Z 0J 


Vfi 

(cu + cu*) 2 + yj 


X(CU + (Of) 




( 16 ) 


One expects that the empirical factor x(cu) approaches unity as cu-*0 and falls off very quickly after 
cu is beyond a certain region. By fitting the experimental data with the expression Eq. (16), one 
is able to obtain an explicit form of the x function. 23 Unfortunately, the spectral density given above 
does not satisfy the detailed balance requirement due to the assumption made on the evenness of 
the x factor. 
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2.4. A modified empirical expression for the spectral density 

One can obtain a modified empirical model by introducing a x factor in the expressions of the 
spectral densities 

f(m) \ (m - (“-“»)• <l7) 


n ~ w) ~ ; £ (a + i r + ~ “ j) - (,8) 

However, instead of assuming that the x factor introduced above is an even function as done 
previously, one imposes the requirement of detailed balance on the spectral density: 

F(-ft)) = e-*“'* r F(cu) 


.1 


= e *“"* 7 -Xp/M 2 ; 


Vfi 


rX(co +w fi ). 


(19) 


*7 (cu + t^r + yj 

Comparing the corresponding resonant terms in expressions (18) and (19), the following equality 
must hold 


x(-co - = e~ H ' a * a f )lkT x{pJ + (o fi ). 

The conclusion is that in order to satisfy the detailed balance requirement, the empirical x factor 
should not be an even function, but 


X(-co) = c~ MkT x(co); 

(20) 

or expressed more symmetrically, 


e? al 2 kT x(.— co) = e~*° >l 2 * T x(co). 

(21) 

Therefore, it is the symmetric factor, x(<o ), defined by 

X(co) = c-^^xico). 

(22) 

and not *(cu) that is an even function; i.e. 

Z(a>) = z(-o>). 

(23) 


In terms of x the modified spectral density can be expressed as 

f (“) - e ‘*'“ r ; | ft*'**" W (M _ ^ 


-c— < 24 > 

where the symmetric line intensities ||/^|| : ( = 11^ II 2 ) are defined by 

IM 2 « (P,P/)'W = g/*-* E ' +El)lu W/Q, (25) 

» 

in which gfi = g,) is the nuclear spin degeneracy factor and Q is the partition function. Similarly, 
using Eqs. (19), (22), (23) and (25), we have 

FI -m) = e--""; I In, F (c ; ■ _ t(w - cjf). ■ (26) 

More clearly, if we define G(a>) by 


G(co) = -Z 


: ft 

(to - (Off + y} 


x ( ct ) - C 0 ji ), 


(27) 
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which is an even function, then the spectral density can be simply expressed as 

F(co) = e»“« r G(ffl), 


F(-oj) = (28) 

Therefore, this modified spectral density does satisfy the detailed balance requirement which can 
be expressed more symmetrically as 

e-^^fa) = e^^Fi-w). (29) 


Finally, the symmetric modified spectral density is 


F(co) + F( — cd) = 2ch(ftco/2kT)G(co) = 5- £ S, 

*«,»o 


ch(hco/2kT) f y fi _ 

ch(hco fi /2kT) |> - o> fi f + y} X(0} a * } 


+ 


yfi 

(w + Oifi) 1 + yj 


X(a +<o fi ) 


— 


(30) 


We note that the modified empirical expression of the spectral density has the same form as that 
derived by Rosenkranz in his theoretical formalism. 1314 

By comparing the above expression with the corresponding expression (16) in which both of the 
empirical factors denoted by z(co) and z(co), respectively, are even functions, one is able to find 
that these two expressions differ from each other by the factor ch(ha)/2kT)/ch(haj fi /2kT). It is 
obvious that in the band center regions where most of the contributions to absorption come from 
nearby strong resonant lines and where a> — a> fi a: 0, the difference between these two formulas can 
be ignored. However, since atmospheric windows are the primary regions in which the empirical 
factors play a crucial role, including or excluding the factor ch (ha>/2kT)/ch {hco fi /2k T) makes 
significant differences. In addition, these differences will increase as the temperature decreases. 
Although the results of continuum absorption obtained from these two formulae can be made to 
coincide in a given frequency region for a specified temperature (for example, in the 10 pm window 
at room temperature for which experimental data are available and on which the empirical fittings 
are partially based), differences will definitely appear when predicting the absorption by interp- 
olation or extrapolation to other temperatures and/or frequency regions; sample numerical tests 
confirm this conclusion. Therefore, the intrinsic weakness of Eq. (16) will limit its reliability for 
atmospheric applications. 


2.5. A more general modified spectral density 

More generally, in order to distinguish the positive-frequency resonance and the negative-fre- 
quency resonance, one has to introduce two z factors, z+ and X- instead of one z factor in the 
empirical model, 19 " 21 


F(<o) = ^- £ 

** atj, >0 


M 2 




(31) 


F(—w) = i Y 

71 Of>0 




j/ V -°» )+ v 4V+ 7 } z - ( - a - +' “'V 


(32) 

Again, imposing the requirement of the detailed balance, F( — co) = e~ tolkT F((o), on the spectral 
density, we must have. 


e*“ /2 * r Z + (~cu) = e - tofl * 7 x_(a>), 

c ka ’ /2Ar X-(-a>) = c- kaatr z+(a>). (33) 


Similarly, we can define 


* + (co) = e- fc ” /a7 x + (co), 

Z-(<o) = e _ *" /J * r x_(<a). (34) 
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We note that neither of these two factors is an even function, rather they are related to each other 
by an inversion operator /, 


= /*-(&>) = Z-( -®), 

Z_(co) = Ix+(a>) = Z+(-co)- 

In terms of them, the modified spectral density can be expressed as 


and 

F(—co) = e~ lm,nkT 
More clearly, if we define 


(35) 


na) - [ ( „ . X'f + y} *♦<" - ^ + c + 1’f + ,} +< *>} (36) 
■* !|1 + ^ + «*<“ - «*>]• (37) 

define 

c ' (ffl) ‘ i l|! + **'“ - 1 "'> ■ + f - ( " + "4 (38) 


which is an even function, then we have 


F((o) = e h °’ niT G'(oj), 

F(-co) = e-^ l ’ 2iT G'(co), (39) 


where the fact that the detailed balance requirement is satisfied is clearly shown. The symmetric 
modified spectral density can accordingly be expressed as 


F(oj) + F(-oj) m 2ch(ha)/2kT)G'((o) 

_ 3 he 1 ^ ch(hco/2kT) [~ y fi . , s 
" 5 * ch(hw fi /2kT) [> -u> fi ) 2 + y} X+{< ° ~ ^ 

+ ^u7/ Am+a ‘> ) } m 

By introducing two x factors, x+(cu) and *_(a)), in the above expression, the asymmetric character 
of the line shapes is explicitly formulated. However, if we do not distinguish between x + and %- 
and replace both of them by the same x, Eq. (40) will be the same as Eq. (30). We note that the 
expression of spectral density in terms of two x factors is not only more general than the expression 
in terms of one x factor, but also has a sound physical basis. 19-21 


3. CONCLUSIONS 

As has been pointed out in Sec. 2.3, the empirical model widely used in transmission codes fails 

to satisfy the detailed balance requirement. Its failure can be traced back to an arbitrary 

assumption that the factor x(a>) introduced to modify the Loren tz line shape is an even function. 

By lifting this assumption and introducing the symmetric empirical function x(co) instead of y(a>), 

one is able to obtain an empirical formalism which satisfies the detailed balance requirement and 

overcomes the intrinsic weakness of the previous formalism. In this formalism, it is ^(co), not 

x(ct>) itself, that is an even function. In addition, in order to derive a more general expression of 

the spectral density, one has to distinguish the positive-resonance terms from the negative-reson- 

ance terms and, consequently, none of the x factors is an even function. Therefore, one has to 

introduce two x factors, one for the positive frequency region and the other for the negative 

frequency region, 3 * * * * * * * * * * * * * * * 19 " 21 in order to take into account the asymmetric character of the line shapes. 

This more general empirical formalism also satisfies the detailed balance requirement and has a 

more sound physical basis. This latter formalism should be used in the further development of 

transmission codes. While this paper has focused on the spectral densities, a parallel study on 

their time domain partners, i.e., the correlation functions has also been performed and will be 

presented elsewhere. 
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